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The exact solution of the spin one-half Falicov-Kimball 
model, with random hopping between the lattice sites, is 
used to explain the anomalous magnetic response of Yb-based 
valence-fluctuating intermetallic compounds. The anoma- 
lous behavior arises from an entropy-driven local-moment- 
nonmagnetic transition of unhybridized Yb ions in these ma- 
terials which can also be used to explain the observed meta- 
magnetism and resistivity anomalies. 



The valence- fluctuating (VF) compounds YbInCu4, 
YbIni_2;Ag2;Cu4 and Ybi_a;Ina;Cu2 exhibit large anoma- 
lies in their thermodynamic spectroscopic ||l|,§j5| 
and transport properties 0,^6 1, around a character- 
istic, sample-dependent temperature Ty (30 — 70K). 
The anomalies are due to the VF, as the L///-edge, 
Mossbauer, X-ray absorption and the thermal-expansion 
data show that the Yb ions fluctuate between a 2+ and 
3-1- state. However, the atomic volume is altered at Ty 
by a small amount and the average valence of the Yb ions 
does not change by more than 4%. The X-ray and 
neutron data do not show any structural changes across 
the VF transition and the neutron Q and the 
NMR data |^ rule out magnetic ordering as a possible ex- 
planation of the anomalous behavior. Here, we focus on 
the anomalous magnetic response of Yb-based VF com- 
pounds, as revealed by the temperature dependence of 
the low-field susceptibility, x(T), and the field- induced 
torque, r(T), and by the measurements of the high-field 
magnetization, M{H), at temperatures below Ty. 

The anomaly in x(T) is seen most clearly in the 
YbInCu4 single-crystal data |3|, as a pronounced asym- 
metric peak just above Ty. The anomalous peak rises 
from a large uniform background and is much steeper on 
the low-temperature than on the high-temperature side. 
Below Ty, x(2^) drops to a minimum and then grows 
slowly towards the lowest temperatures, where a weaker 
Curie- like upturn is often seen. For T 3> TV, x{T) de- 
cays in a Curie- Weiss fashion but the magnitude of the 
susceptibility at 300 K is about the same as at the low- 
temperature minimum. Surprisingly, the magnitude, and 
even the sign, of the Curie- Weiss temperature appears to 
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be sample dependent. Such behavior is quite different 
from what one finds in a typical metallic heavy-fermion 
sample, like YbAgCu4 ]lO|,p[, where the low-temperature 
response is much enhanced over the room-temperature 
one and x{T) exhibits a shallow and symmetric maxi- 
mum. 

The torque induced by the field (measured for the sin- 
gle crystals used in Ref. ) shows that the magnetic sus- 
ceptibility is also anisotropic j^, which is quite unusual 
for compounds of cubic symmetry. The torque, F(r), 
is characterized by an asymmetric peak, which is clearly 
related to the anomalous peak in x{T)- Below Ty, r(T) 
drops by several orders of magnitude and then, like x(T), 
it shows an upturn as the temperature is further reduced. 
The decay of F(T) above Ty is very rapid and cannot be 
described by a Curie- Weiss law. The F-anomaly, like the 
X-anomaly, is sharper and more asymmetric for the sam- 
ples with lower Ty's. 

The anomaly in the high-field magnetization of 
YbInCu4 and YbIni_a;Aga;Cu4 is characterized by a 
sudden increase of the slope and the saturation of M [H) 
at about Hy — 30 — 50 T. The magnetostriction data ||] 
indicate that the metamagnetism relates to valence fluc- 
tuations of Yb ions but the field- induced change in the av- 
erage f-count, {uf) is even smaller than the temperature- 
induced one. At low temperatures, the Zeeman energy 
at Hy is comparable to the thermal energy at Ty. At el- 
evated temperatures, the metamagnetic transition is ab- 
sent [|. The magnetic anomalies at Ty and Hy, and 
the characteristic energies ksTy and fisHy are sample 
dependent while the low-field response, for T <Ty, 
is quite robust with respect to doping and annealing. 

The simultaneous analysis of the susceptibility and the 
torque data shows Q that the magnetic response is due 
to two physically distinct components. The first one, 
XVF, is large (on an absolute scale), isotropic, and sam- 
ple independent. The second one, xiocai, is vanishingly 
small below Ty and smaller than xvf at high tempera- 
tures, but dominates x(T) and r(T) near Ty. xiocai is 
anisotropic and strongly sample dependent. In addition, 
the field and sample dependence of the M{H) data also 
supports the picture of two types of magnetic excitations. 

To explain these features we recall that the Yb com- 
pounds mentioned above crystallize in C15&-(AuBe5)- 
type structure with In and Yb ions at the so-called 
4a and 4c sites, respectively. Because of the short 4a- 
4c distance, the f-states of the Yb ions hybridize with 
the In ions and give rise to an enhanced background sus- 
ceptibility, Xvf -I which dominates x{T) much above and 



1 



below Ty- For T < Ty, these hybridized states also yield 
the sample-independent low-field magnetization. On the 
other hand, the disordered CIS-type structure, in which 
the 4a and 4c sites are randomly occupied by Yb and In 
ions, is quite close to the C15b-type structure, and is not 
ruled out by the X-ray and neutron data (especially 
for a small number of disordered Yb ions occupying the 
In 4a sites). Thus, we assume that a small number (N) 
of Yb ions in YbInCu4 occupy the In 4a sites and that 
their 4f-states remain unhybridized. The hybridization 
between these ill-placed 4a- Yb ions and the regular 4a-In 
ions is absent because the 4a-4a distance is nearly twice 
the 4a-4b distance. The magnetization, of the 4a- Yb ions 
in the f^^ configuration is angle-dependent, because the 
g-factor of the unhybridized /-holes becomes anisotropic 
in the distorted crystal environment. The external field 
induces the torque, which relates to the local susceptibil- 
ity as, r ~ {xiocaiiN, T)) H^, where Ag^ ^ gl ~ g^ 
is the 51-factor anisotropy between the x and y direction, 
and (• • •) denotes the angular average. As show below, 
the unhybridized Yb ions crossover at Ty from the mag- 
netic 3-|- to the nonmagnetic 2+ configuration, so that 
Xiocai and r vanish. At low temperatures and high fields 
the unhybridized states lead to the metamagnetic transi- 
tion with a small energy difference between the low-field 
and the high field states. Note, (uf) does not change 
much either at Hy or at Ty, because the VF transition 
involves only N 4a-impurities out of a total of Ni Yb 
ions, i.e. N/Ni < 1. 

Our model for the anomalous magnetic response of 
these materials begins with this assumption that the f- 
electrons of the disordered Yb ions are not hybridized 
with the rest of the lattice. These ions form a random 
lattice. The conduction electrons can move, via nearest- 
neighbor hopping on the regular lattice, between any of 
these random lattice sites; i.e. the system can be modeled 
by localized f-states at energy E and delocalized d-states, 
with a random hopping between each of the sites of the 
lattice. We assume that the / and d states have a com- 
mon chemical potential (/i) and interact by a Coulomb 
repulsion (C7). In addition, both the d and / particles 
carry a spin label a and the d-level can accommodate 2 
electrons (or holes) of the opposite spin. The occupancy 
of the /-level is restricted to n/ < 1 because of the large 
Coulomb repulsion {Uff ^ 00) of the /-particles of op- 
posite spins. To discuss the Yb-based VF compounds 
we use the hole-picture, in which E < fi and the to- 
tal number of holes at the ill-placed sites is restricted 
to n/j = -f < 3. In the electron-picture, one has 
E > fi and restricts the total number of electrons to 
Ue = n"^ + n'j < 3. The magnetic field h couples to the 
/ and d states but with different 5-factors {g — 4.5 for 
the /-holes) . This picture is described by a spin one- half 
Falicov-Kimball (FK) model 
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where the notation of Ref. is used and i and j denote 
the sublattice sites where the ill-placed Yb ions are lo- 
cated. The hopping matrix elements are chosen from 
a random distribution in such a fashion that as A'^ ^ 00 
the noninteracting density of states becomes Gaussian 
p(e) = t* exp[-e2/i*2]/0F ||l3|. The only role of the reg- 
ular lattice, with A'^; sites, is to provide an infinite-range 
(random) connectivity between the sublattice sites. 

In the hmit A^ ^ cxo Q, and N/Ni < 1, our choice 
for tij maps this problem onto the infinite-dimensional 
(local) one, which allows the magnetic susceptibility to 
be evaluated exactly using the methods developed by 
Brandt and Mielsch |12|. The only differences we have is 
that we fix the total electron concentration (by adjust- 
ing /i) rather than fixing just the /-electron concentra- 
tion. The Green's functions and susceptibilities are de- 
termined numerically by solving the self-consistent equa- 
tions of Ref. ||l2|. The properties of the model depend 
on n/i, i? — /i, and U] the value of the effective ma- 
trix clement for random hopping (t*) defines the energy 
scale. The uniform f-f spin-spin susceptibility satisfies 
X^SDwiT) = n)[T)/2T. 

We present the susceptibility results of our numeri- 
cal solutions for a variety of different cases in Figure 1. 
Each figure corresponds to a different value of n/i, with 
ten values of E/t* ranging from to —4.5 in steps of 
—0.5. Fig. 1(a) is a typical high-density result. There 
are 2.5 holes per impurity site, implying 7^ at all 
temperatures, leading to a Curie contribution to the sus- 
ceptibility at small T. In this regime, the results are 
rather insensitive to U [U — IQt* here). Fig. 1(b) is the 
Uh = 2.1 case. We chose U — t* here, and the results de- 
pend sensitively on E, showing a downturn at moderate 
T, before the Curie-law divergence sets in for T — 0. As 
n/i is decreased to 2 and beyond, it is no longer neces- 
sary for there to be any / holes remaining at T = 0, and 
X-^-^ can vanish in that limit. In Fig. 1(c) we plot X'^^ 
for Uh — 2 and U — 2t* , showing an asymmetric peak 
that increases in size and sharpens as Ty decreases, and 
which still has a weaker Curie upturn at the lowest tem- 
peratures, because ^ as T ^ 0. Here, the results 
depend quite sensitively on E and U , with the peak low- 
ering in magnitude and broadening as U increases, and 
the low-temperature upturn becoming more prominent. 
The case Uh = 1.9 is shown in Fig. 1(d) for U — lOt*. 
Once again we have a sensitive dependence on the param- 
eters, but the low-temperature upturn has disappeared. 
In some cases X'^^ has a double-hump structure. Finally 
the low-density regime is plotted in Fig. 1(e) {uh = 0.5, 
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U = lot*). Here, as in the high-density hmit, the results 
are insensitive to U . 

The model also exhibits a metamagnetic transition, for 
large enough magnetic fields. This is shown in Figure 2 
for one case: nn = 2, E = —t*, U = 2t*, and eight 
temperatures from 0.05t* up to 6.4i*. The numerical 
analysis shows that at T = 0, the Zeeman energy at the 
transition is of the order of fcsTy . Hy decreases at higher 
temperatures but the the transition becomes smoother. 
The metamagnetism disappears as T is increased beyond 
Ty, because the /-holes are already occupied in the zero- 
field limit. 

The exact solution reproduces well the overall behav- 
ior of the experimental data. Figs.(lc)-(ld) capture most 
of the features shown by x{T) in Refs. 0-ii,0,§, 

such 

as an asymmetric peak at Ty that increases in magni- 
tude and sharpens as Ty 0, while Fig. (2) explains 
M{T) reported in with a metamagnetic transition 

that smoothes out and then disappears as T is increased. 
Various samples will have different numbers of Yb im- 
purities and will require different values of the coupling 
constants. Our analysis shows that the lower the transi- 
tion temperature, the more pronounced and steeper the 
anomaly. For large fields, unhybridized Yb ions crossover 
to a magnetic configuration at much lower temperatures, 
than in the absence of the field. 

Although this is a complicated many-body system, in 
which the behavior depends strongly on the total number 
of particles at the impurity sites and the double occu- 
pancy of 4/-states is prohibited, the anomalous behavior 
can be understood in terms of simple thermodynamic 
considerations. For E < fi and rih < 2, the nonmag- 
netic state of an Yb impurity is energetically more favor- 
able than the magnetic state, so the ground state corre- 
sponds to the 4/^"* configuration (no /-holes). However, 
at high T, the large entropy of the /-hole spins favors 
the magnetic state. Thus, the high-temperature response 
of an ensemble of Yb impurities approaches the limit 
XiocaiiN,T) ~ NxoiT), where xo(T) ~ {giiB? /T is the 
local susceptibility of a single Yb ion with one unhy- 
bridized magnetic hole. The g-factor anisotropy of that 
hole gives rise to the field-induced torque which depends 
on temperature via xiocai{T). The anomalous part of the 
susceptibility and the torque grow rapidly as the temper- 
ature is reduced until, close to Ty, the entropy gain is in- 
sufficient to compensate the energy loss and the / holes 
are filled to form a nonmagnetic state. Thus, xiocai drops 
significantly, which is the origin of the sharp asymmetric 
peak in the r(T) and x(T') data. This anomalous peak 
conceals the broad maximum which is usually found in 
the susceptibility of VF compounds (from the electrons 
on the regular lattice). For rih > 2, however, the /-hole 
cannot be completely filled at T = and the residual /- 
spin leads to a low-temperaturc upturn of x{T), which is 
also often seen in the data ||l|-^,^j^J^ (this theoretically 
predicted Curie tail will add to the Curie tail present in 



all samples due to impurities). The associated Curie con- 
stant is very small for n/j ~ 2 but increases rapidly as the 
density of holes increases. In YbInCu4, the proximity of 
the f-level to gives rise both to the FK-transition of 
unhybridized Yb ions and the VF behavior of hybridized 
Yb ions. Had E — ^ been much larger than t* , then the 
unhybridized f-levels would not change from 3+ to 2+ at 
low temperatures. 

The transport properties are also explained by this 
scheme: the Kondo scattering for T > Ty ^ Tk is at the 
spin-disorder limit, while for T < Ty the free spins vanish 
and the Kondo scattering is absent. Hence, the resistivity 
jumps at Ty from its small low-temperature value to a 
high-temperature spin-disordered value [^. The switch- 
ing on and off of the Kondo scattering also leads to large 
anomalies in the elastic constants through the electron- 
phonon coupling at the regular lattice sites. 

In summary, the response of YbInCu4 and other Yb- 
based VF compounds is analyzed in terms of two inde- 
pendent magnetic excitations. The magnetic anomalies 
are attributed to an entropy-driven transition of disor- 
dered Yb and explained by the Falicov-Kimball model 
with random hopping. Although the disordered Yb ions 
undergo at Ty a large valence change, (n/) does not 
change much since N/Ni ^ 1. Thus we account for the 
correlation between x(^) r(T), explain the meta- 
magnetic transition, and reconcile the large changes in 
magnetic, transport, and elastic properties with small 
changes in the /-count and the average atomic volume. 
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FIG. 1. Magnetic susceptibility of the /-holes in the Fal- 
icov-Kimball model. The different curves correspond to 
E/t* = 0.0, —0.5, . . . , —4.5 (in general E decreases from top 
to bottom in these figures). The different figures are: (a) 

n** = 2.5 and U = lOt* ; (b) n'^ = 2.1 and U = t*; (c) n'' = 2.0 
and U = 2t*; (d) n'' = 1.9 and U = lOt*; and (e) n'' = 0.5 
and U = Wt*. 
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Freericks and Zlatic, Phys. Rev. Lett., Figure 1(b) 
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Freericks and Zlatic, Plnys. Rev. Lett., Figure 1(c) 



FIG. 2. Total magnetization of the Falicov-Kimball model 
as a function of magnetic field. The different curves corre- 
spond to different temperatures: 
T/t* = 0.05,0.1,0.2, ...,3.2,6.4 (the temperature increases 
from top to bottom in the large- /i range). The parameters 
are = 2, E = -t* , and U = 2t* . 



Q 




Freericks and Zlatic, Phys. Rev. Lett., Figure 1(a) 
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Freericks and Zlatic, Phys. Rev. Lett., Figure 1 (d) 
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Freericks and Zlatic, Phys. Rev. Lett., Figure 1 (e) 



4 




h/t* 

Freericks and Zlatic, Phys. Rev. Lett., Figure 2 



